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Abstract. The natural spherical projection associated with the Hierarchical Equal Area and isoLatitude Pixelisation, HEALPix, 
is described and shown to be one of an infinite class not previously documented in the cartographic literature. Projection 
equations are derived for the class in general and it is shown that the HEALPix projection suggests a simple method (a) 
of storing, and (b) visualising data sampled on the grid of the HEALPix pixelisation, and also suggests an extension of the 
pixelisation that is better suited for these purposes. Potentially useful properties of other members of the class are described. 
Finally, the formalism is defined for representing any member of the class in the FITS data format. 

Key words, astronomical data bases: miscellaneous - cosmic microwave background - cosmology: observations - methods: 
data analysis, statistical - techniques: image processing 



1. Introduction 

The Hierarchical Equal Area and isoLatitude Pixelisation, 
HEALPix (Gorski et al. 2004) offers a scheme for distributing 
Y1N^(N € N) points as uniformly as possible over the surface 
of the unit sphere subject to the constraint that the points lie 
on a relatively small number (4A^ - 1) of parallels of latitude 
and are equally spaced in longitude on each of these parallels. 
These properties were chosen to optimise spherical harmonic 
analysis and other computations performed on the sphere. 

In fact, HEALPix goes further than simply defining a dis- 
tribution of points; it also specifies the boundary between ad- 
jacent points and does so in such a way that each occupies the 
same area. Thus HEALPix is described as an equal area pix- 
elisation. Pixels at the same absolute value of the latitude have 
the same shape, though pixel shape differs between latitudes. 
The boundaries for N - I define the 12 base- resolution pix- 
els and higher-order pixelisations are defined by their regular 
subdivision. Note, however, that although they are four-sided, 
HEALPix pixels are not spherical quadrilaterals because their 
edges are not great circle arcs. 

HEALPix was originally described purely with reference 
to the sphere, the data itself being stored as a one-dimensional 
array in a FITS binary table (Cotton et al. 1995) with either 
ring or nested organization, the former being suited for spher- 
ical harmonic analysis and the latter for nearest-neighbour 
searches. For visualisation purposes the software that imple- 
ments HEALPix offers a choice of four conventional projection 
types onto which HEALPix data may be regridded. 
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Fig. 1. The HEALPix class of projections for H = 1,2,3 (top), and 
the nominative case with H = A (bottom) at twice the linear scale and 
with the top left-hand comer of the graticule "cut away" to reveal the 
underlying cylindrical equal-area projection in the equatorial region. 
Facets are shown as dashed diamonds. 



However the existence of the analytical mapping of each 
of the twelve base-resolution pixels onto a [0, 1] x [0, 1] unit 
square was noted, and in fact the HEALPix software uses this. 
Roukema & Lew (2004) have recently provided a mathemati- 
cal derivation of the equations and present a diagram showing 
a projection of the whole sphere (hereinafter the HEALPix pro- 
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jection) in which the base-resolution pixels, and consequently 
the pixels of all higher-order pixeUsations, are projected as dia- 
monds (i.e. squares rotated by 45°). These equations may read- 
ily be synthesised into those of an equal area projection of the 
whole sphere. 

The HEALPix projection does not appear to have been doc- 
umented previously in cartography texts and could not be lo- 
cated in a web search; in particular, it is absent from Snyder's 
(1993) review of the history of cartography, and also from 
Pearson (1990). Ilhnois State University's MicroCAM web- 
site presents a catalogue of 320 map projections produced 
by a member of the International Cartographic Association's 
Connmission on Map Projections (Anderson, 2003). None bear 
even a superficial resemblence; the equal-area quad-cube may 
be disected and rearranged to produce something with a simi- 
lar boundary but it is a distinctly differentprojection. The stated 
intention of this website is to present as complete a collection 
as possible of historical, published map projections. 

This work will show that the HEALPix projection is one of 
the more important members of an infinite class of projections 
parameterised by // € N and will derive the projection equa- 
tions for the class. In particular, the HEALPix projection (i.e. 
with H - A) suggests a simple way of storing HEALPix data on 
a two-dimensional square grid as used in conventional imaging 
and mapping, and also suggests an extension of the HEALPix 
pixelisation that is better suited to this. The HEALPix projec- 
tions with H = 3, and 6 are also shown to be special, and their 
properties will be discussed. 

The related issues of representing celestial coordinates in 
the HEALPix projection are also considered in relation to im- 
age data storage in FITS (Hanisch et al. 2001). 



2.1.1. Transition latitude, 6x 

In deriving the projection equations, note firstly that for any 
H the total area occupied by the half-facets in the north polar 
region is always 1/6 of the total area. Since the projections are 
equi-areal, we equate the area of a spherical cap on the unit 
sphere, A = 2n{\ - smO), with the corresponding fraction of 
the total area, 47r/6, to obtain the transition latitude, 6^, which 
is independent of H: 



Ox = sin"'(2/3) « 41°8103. 



2.1.2. Equatorial region 



(1) 



The equatorial region, where \6\ < Ox, is clearly a cylindrical 
equal-area projection, i.e. {x, y) = (0, a sin (9), where a is a con- 
stant determined by the requirement that be projected at the 
vertex of a square, i.e. - n/H = a sin 6^, whence 



X - (p, 

3n . 
u = — sm 6. 
^ 2H 

Because dy/d(f> 
Jacobian reduces to 



J{<f>,e) = 



1 dx dy 
co&ed^W 



(2) 
(3) 

for the HEALPix projections the 



(4) 



This gives the ratio of an infinitesimal area in the plane of pro- 
jection to the corresponding area on the sphere. In the equato- 
rial regions it is 3n/2H, a constant, indicative of an equi-areal 
projection. Note that the Jacobian is inversely proportional to 
H, but the graticules in the top part of Fig. 1 were set to the 
same areal scale by scaling x and y by ^fH. 



2. The HEALPix projections 

Firstly in this section we derive the projection equations for the 
infinite class of HEALPix projections. This leads naturally to 
the subject of mapping HEALPix data on a regular grid and 
extensions of the HEALPix pixeUsation that arise. 

2.1. HEALPix derivation 

Figure 1 shows the first four members (H - 1 , . . . , 4) of the 
HEALPix projections. They may be described as interrupted, 
equal area, pseudo-cylindrical projections whose defining char- 
acteristics are 

1. They are equi-areal; regions with equal areas on the sphere 
have equal areas in the plane of projection. 

2. Parallels of latitude are projected as horizontal straight fines 
(interrupted in the polar regions) whence dyjdcj) = 0. 

3. Parallels are uniformly divided (apart from interruptions). 

4. The interruptions are defined by stacking equal-area dia- 
monds (hereinafter facets) as shown in Fig. 1. The facet 
that straddles ±180° is split into halves in the graticule. 



2.1.3. Polar regions 

In the polar regions the area north of 6 (> 6y) on the unit sphere 
is 27r(l - sin0) and, noting that the pole is projected at y - 
2njH, in the plane of projection it is H(2n/H - y)^. Equating 
the ratio of these to the value of the Jacobian obtained for the 
equatorial region and solving we obtain 

y = ±^i2-(r), (5) 

where the negative sign is taken for the south polar region, and 

cr = V3(l - Isin^l) (6) 

is the ratio of the distance of the pole from the parallel of 6 to 
that of the pole from the parallel of Ox ■ 

The equation for x may be obtained readily by integrating 
Eq. (4) with dy/dO from Eqs. (5) & (6) to produce x - o-(p + C, 
where C is the constant of integration, thus indicating that the 
parallels are uniformly divided. It is instructive also to consider 
a geometrical argument; the area of any triangle in the ix,y) 
plane with its apex at the pole and base along a given parallel 
of latitude depends only on the change in x between its base 
vertices and not on their location. Since the projection is equi- 
areal, x must therefore vary linearly with (p. 
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Applying the interruptions to the parallels (which in fact 
could be omitted or done in other ways to produce different 
projection types) we have 



X = 0c + (0 - 0c) cr, 
where 

+ 7t)H 



00 = ^ 2 



In 



+ 1 



(7) 



(8) 



f dx dy 
f dx dy 



, 37r 

1, — cos 
2H 



= cr. 



37T cos 6 
2H o- 



. . . equatorial 
. . . polar 



(9) 
(10) 



which shows that in the polar regions x is scaled directly, and y 
is scaled inversely by cr{ff) in order for the Jacobian to maintain 
constancy. 

To get some idea of the relative degree of distortion be- 
tween members of the class, consider first from Eqs. (3) & (5) 
that y scales asl/H for any 6, while x is independent of H ex- 
cept for defining the interruptions. Hence the relative spacing 
of parallels between the equator and poles is independent of H, 
as is evident in Fig. 1, and the distortion is determined solely 
by the relative y : x scaUng. 



is the native longitude in the middle of a polar facet and [mJ , 
the floor function, gives the largest integer < u. 

2.1.4. Properties 

The most important feature of the HEALPix projections, in- 
deed the underlying rationale for the HEALPix pixelisation, 
is that they are equi-areal with squared boundaries. Thus they 
may be completely inscribed by diamonds of equal area, the 
minimum number of which is 3H (the facets). Each facet is 
subject to further subdivision by A^^ smaller equal-sized di- 
amonds that are identified as pixels; their centre positions in 
((/>, 6) may be computed readily for any {H, N) from the inverse 
of the above projection equations as are cited in Sect. 2.4. As 
explained by Gorski et al. (2004), it is significant for spheri- 
cal harmonic analysis that the pixel centres he on a relatively 
small number of parallels of latitude, and that the facets may 
be subdivided in a hierarchical way. 

Of course a pixehsation may be constructed similarly from 
a cyhndrical equal-area projection, but the HEALPix projec- 
tions are much less distorted in the polar regions than any 
such projection. Consequently the HEALPix pixels are much 
truer in shape when projected onto the sphere and their centres 
are much more uniformly distributed. As shown by the dashed 
lines in the upper- left corner of Fig. 1 , the equivalent portion of 
the underlying cylindrical projection, being severely squashed 
at the pole, is stretched upwards to twice its height and brought 
to a point; the pole itself is thereby projected as H points rather 
than a Une. However, this is gained at the cost of introducing 
H — \ interruptions which should properly be considered as 
extreme distortions, though of httle consequence for the pixeli- 
sation. 

Evaluating the partial derivatives we find 




HEALPix image grid 



Fig. 2. The HEALPix pixehsation for Af = 6 on the HEALPix projec- 
tion for = 4 projected with a 45° rotation onto the mapping grid 
showing the twelve facets with a convenient numbering scheme. The 
graticule of the HEALPix projection is shown in the seven facets adja- 
cent to (0, (?) = (0, 0), and those at lower left show the pixel boundaries 
for Af = 6 as defined by the HEALPix pixehsation. 

A spherical projection is conformal or orthomorphic (true 
shape) at points where the meridians and parallels are orthog- 
onal and equi-scaled. The general equations of the cylindrical 
equal area projection expressed in terms of the conformal or 



standard latitude, 6o, are {x,y) 



, sin 61 cos 6o) (e.g. see 



Sect. 5.2.2 of Calabretta & Greisen 2002), whence from Eq. (3) 



-1 l^H 

9o = cos -i/ 

V 3n 



(11) 



For // = 1, 2, 3, 4 this is (62?57, 49°35, 37°07, 22?88); the first 
two of these exceed Ox and hence are inadmissible, and 6o is 
undefined for higher values of H. Since the latitude that halves 
the area of the equatorial region is sin"'(l/3) = 19°A7, inde- 
pendent of H, this suggests that the projection with H = 4 h 
the least distorted in the equatorial regions. 

Looking at it another way, the requirement for equiscaling 
in X and y where the meridians and parallels are orthogonal, i.e. 
everywhere in the equatorial region, and along the centreline in 
the polar half-facets, is 



1 dx 



dy 
do' 



cos 6 d(f> 

Substituting Eqs. (9) and (10) gives 



— cos^ e 

2 



Ho = -(1 -Hlsin^l) 



. . equatorial 
. polar, centreline 



(12) 

(13) 

(14) 



which gives i/o = (4.7,4.4,3.5,2.6,2.7,2.9,3.1,3.1) for 6» = 
(0, 15, 30, ffx, 45, 60, 75, 90). Thus = 4 is a good all-over 
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Fig. 3. HEALPix double-pixelisation for N = 3 on the HEALPix pro- 
jection with H = 4. Filled circles define the regular grid, with inter- 
polated pixels shown as open circles to the left of these. Pixel bound- 
aries are shown in the middle three facets, those of the two additional 
polar pixels contain a contribution from each of the four adjacent po- 
lar facets. Note the difl^erence between the pixel locations in this fig- 
ure compared to the Af = 6 pixeUsation in Fig. 2. Also illustrated is 
the boundary between the faces of the pseudo-quadcube layout of the 
HEALPix projection (grey). 

compromise but for \9\ > 30°, the latitude that halves the area 
of the hemisphere, H = 3 would appear to be better on this 
basis. 

The nature of the projective distortion in the region where 
meridians and parallels are not orthogonal is more complicated. 
In the polar regions the projection of the facets onto the sphere 
(i.e. the base-resolution pixels) meet at the pole at 360°///. For 
H = 4 this is 90° which accords with the angle in the plane 
of all HEALPix projections. Thus it might seem that H = 4 
should be least distorted in the neighbourhood of the pole. 
However, this argument is specious; on the sphere the angle 
between meridians and parallels along the edges of the polar 
half-facets is always 90°, while in the plane of projection it is 
always 45°. 

These coimnents on distortion will be qualified in Sect. 2.3 
for members of the class for which additional j/-scaling is ap- 
plied. 

2.2. The HEALPix grid 

The base-resolution pixels of the HEALPix pixelisation are 
projected as diamonds (squares rotated by 45°) on the 
HEALPix projection with the consequence that the pixel lo- 
cations fall on a grid with diamond- symmetry. 

However, Fig. 2 shows that the diamond grid may be con- 
verted to the common square grid used in imaging via a triv- 
ial 45° rotation. The resulting image plane is slightly less than 
half-filled (48%) but this is comparable to the figure of 50% for 
quad-cube projections (O'Neill & Laubscher 1976) which are 
commonly used in the same type of application as HEALPix. 
Moreover, being composed of square facets Uke the quad- 
cubes, the HEALPix projection also admits the possibility of 
dissection and storage on a third image axis, such as is imple- 




Fig. 4. HEALPix projection for H = 6 scaled in !/ by V3 whereby the 

square facets become pairs of equilateral triangles. These may then be 
subdivided in a hierarchical way; base-level pixels are shown to the 
left and right (black circles), the first level of subdivision is at mid-top 
(grey circles), the next at mid-bottom (open circles). 



mented for the quad-cubes via the CUBEFACE keyword in FTTS 
(Calabretta & Greisen (2002). In this context Fig. 3 shows how 
the facets may be repartitioned into a configuration that resem- 
bles that of the quad-cube faces. However, this resemblence is 
purely superficial because the "cubeface" edges do not match 
those of a quad-cube projection on the sphere. 

Facet number 4 which straddles (p = ±180° may be treated 
in a number of ways; it may be left split, or the halves may be 
reconnected in either the lower-left or upper-right comer, or it 
may be replicated in both. 

2.2.1 . HEALPix double-pixelisation 

The main drawback with the above technique for storing 
HEALPix imaging data is that the image is presented at an 
unusual orientation. However, this may be solved via a sim- 
ple extension to the HEALPix pixelisation. Figure 3 shows the 
HEALPix grid with a pixel interposed between every pair of 
pixels along the parallels of latitude and additional pixels added 
at the two poles. The total number of pixels in the pixelisation is 
thereby increased from 12A^^ to 24 A^^ + 2 without affecting the 
special properties described by Gorski et al. (2004), although 
requiring a slightly different method of forming the hierarchy 
and indexing it. Pixels that fall along the fines where the polar 
half-facets meet are distorted in such a way as to "zip" the two 
edges together but they still have equal area. 

2.3. Other pixelisations 

Consider dividing the 360° of circumpolar latitude into integral 
subdivisions. Of the possible ways of doing this (1 x 360°, 2 x 
180°, 3 X 120°,4 X 90°,5 X 72°,6 x 60°, . . .) only the divisions 
into 3, 4, and 6 correspond to regular polyhedra. The division 
into 4 X 90° corresponds to the famiUar case of HEALPix with 
H = 4 with diamonds tesselated by diamonds. 

2.3.1. Triangular - 7/ = 6 

However, the division into 6 X 60° suggests a different type 
of pixeUsation in which equilateral triangles are tesselated 
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Fig. 5. HEALPix projection for /:/ = 3 scaled in !/ by 1/ V3 whereupon 
it becomes three consecutive hexagons. 

by equilateral triangles. This pixelisation may be defined by 
rescaling the HEALPix projection with = 6 by V3 in j/ 
so that the half-facets become equilateral triangles. Such a lin- 
ear scaling does not affect the projection's equal area property. 
What were previously half-facets may now be identified with 
36 new, triangular base-resolution pixels that may be subdi- 
vided in a hierarchical way as for HEALPix, as depicted in 
Fig. 4. It is interesting to note that this subdivision is natu- 
rally hierarchical - the number of pixels varies exponentially as 
36 X 4^^"' where A' is the hierarchy level. In the H = 4 pixelisa- 
tion the exponential hierarchy must be engineered by doubling 
N. 

The conformal latitude computed for H - 6 with this 
extra (/-scaling is 9o = 30°98, indicating that the projec- 
tion becomes conformal at the latitude that bisects the hemi- 
sphere by area. Applying Eqs. (13) and (14) with the extra 
scaling gives ^/3H, = (8.2,7.6,6.1,4.5,4.6,5.1,5.3,5.4) for 
6 = (0, 15, 30, ^x, 45, 60, 75, 90), again indicating less distor- 
tion in the polar regions than // = 4. It also does better in 
the polar zone away from the centreline because the 60° angle 
along the edge of the polar facets more closely matches the true 
angle of 90° on the sphere. Overall, this pixelisation performs 
adequately at low latitudes and does better than the H = 4 pix- 
elisations at mid to high latitudes. 

This rescaling of the H - 6 projection is reminiscent of 
Tegmark's (1996) icosahedral projection composed of 20 equi- 
lateral triangles; the problem of indexing the subdivisions of its 
triangular facets was solved in the implementation of the cor- 
responding pixelisation. In the present context the iso-latitude 
property is still present but modified somewhat from the dia- 
mond pixelisation ofH = 4. However, if the pixel centres are 
moved up or down from the centroid by of the height of 
the equilateral triangles to the point half-way between the base 
and apex then they fall onto a regular grid sampled more fre- 
quently in X than y. This provides some of the same benefits as 
the H = 4 double-pixelisation. 

However, although the displacement is small, there is a pos- 
sibility that it will introduce statistical biases so the full con- 
sequences should be investigated for a particular application. 
These potential biases may be minimised by making the pixel 
size sufficiently small, and the fact that the bias between ad- 
jacent pairs of pixels is in opposite senses will tend to cancel 
them over a region encompassing a sufficient number of pixels. 
It should also be remembered that although the pixel locations 
appear to be at the centre of the pixel boundary in the projec- 
tion of the diamond, square, and triangular pixelisations this 
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is very much an artifact of the distortions inherent in the pro- 
jection. Because the ^-coordinate varies non-linearly with 6, 
on the sphere they are actually biased to one side of the pixel. 
Hence some degree of bias is unavoidable. 

2.3.2. Hexagonal -H^3 

The division into 3 X 120° suggests hexagonal base-resolution 
pixels. Although the familiar "honeycomb" structure shows 
that it is possible to tile the plane with hexagons, neverthe- 
less there is no bounded tesselation of hexagons by hexagons. 
That is, a hexagonal region may not be cut out of a honey- 
comb tesselation without cutting the individual elements. Thus 
it may seem surprising that a hexagonal pixelisation can be con- 
structed from the HEALPix projection for H -T> with y scaled 
by 1 / V3. The boundary of this projection, as seen in Fig. 5, is 
reduced to that of three sequential hexagons. This boundary is 
then used conceptually as a "pie-cutter" on a honeycomb tesse- 
lation of the right scale. Pixels that are cut can be made whole 
again by borrowing from adjacent facets, much as the square 
pixelisation in Fig. 3 does. 

Rescaling Eqs. (13) and (14) gives i/o/VI = (2.7,2.5, 
2.0, 1.5, 1.5, 1.7, 1.8, 1.8) for 61 = (0, 15,30, 6lx, 45, 60, 75, 90). 
Thus the rescaled H - 3 projection does not achieve confor- 
mality at any latitude, it does well close to the equator, but 
degrades at mid-latitudes. In the polar regions the 30° angle 
between meridians and parallels along the edge of the facets 
is further from the ideal of 90° than the 45° for the unsealed 
projections. 

2.4. HPX; HEALPix in FITS 

In this section the HEALPix projections are described in the 
same terms as the projections defined in Calabretta & Greisen 
(2002). 

HEALPix projections will be denoted in FITS with algo- 
rithm code HPX in the CTYPE ia keywords for the celestial axes. 
As data storage has become much less of an issue in recent 
years we do not consider it necessary to create an analogue of 
the CUBEFACE keyword to cover HPX. However, if HEALPix 
data is repackaged into the pseudo-quadcube layout shown 
in Fig. 3 the CUBEFACE storage mechanism is applicable for 
H = 4 and will be treated properly by wcslib (Calabretta, 
1995). 

Since the HEALPix projections are constructed with the 
origin of the native coordinate system at the reference point, 
we set 

(0O,eo)HEALPix = (O,O). (15) 

The projection equations and their inverses, re-expressed in the 
form required by FITS, are now summarised formally. 
In the equatorial zone where | sin 1 < 2/3: 

X = 0, (16) 
270° 

?/=— sine, (17) 
in the polar zones, where | sin 6* | > 2/3: 

X = (f>c + ((p - 4>c) 0-, (18) 
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y 



180°^ 
.^(2-0. 



(19) 



where the positive sign on y is taken for 6 > 0, negative other- 
wise, and 



cr = V3(l-|sin0|), 



<f>c = 



180° 



{(/> + m°)H 



360° 



+ 1 - 180°. 



(20) 
(21) 



These equations are readily invertible. In the equatorial 
zone where I !/ 1 < 180° /H: 



= sm 



)■ 



\270° 

in the polar zones, where \ y\> 180° /H: 
4> = + (x - Xc)/cr, 

e = ± sin"' 1 1 



(22) 
(23) 

(24) 
(25) 



where the positive sign on is taken for y > 0, negative other 
wise, and 

\yH\ 



o- = 2- 



180°' 



180° 



{x + no°)H 



360° 



+ 1-180°. 



(26) 

(27) 



where Xc is the value of x in the middle of a polar facet, as for 

<Pc- 

FITS keyword PV/.la attached to latitude coordinate i will 
be used to specify H with default value 4. 

HPX has been implemented in version 3.7 of wcslib which 
is distributed under a GNU Library PubUc License (GLPL). 
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